Wave mixing of hybrid Bogoliubov modes in a Bose-Einstein condensate 
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Mode-mixing of coherent excitations of a trapped Bose-Einstein condensate is modelled using 
the Bogoliubov approximation. Hybridization of the modes of the breather (I — 0) and surface 
(I = 4) states leads to the formation of a Bogoliubov dark state. Calculations are presented for 
second-harmonic generation between the two lowest-lying even-parity m = modes in an oblate 
spheroidal trap. Two hybrid modes are strongly excited near second-harmonic resonance, and the 
coupling strength, and hence conversion rate, to the breather mode is half that given by an equivalent 
hydrodynamic estimate. 
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Observations of coherence of matter wave fields such 
as four- wave mixing Q , squeezing Q and harmonic gen- 
eration 3] , in atomic Bose-Einstein condensates have re- 
cently been reported. The well-defined phase of the con- 
densate means that small amplitude quasiparticle exci- 
tations can be produced that combine, through elastic 
coherent (phase preserving) atomic collisions, to produce 
frequency-mixed modes. In this Letter we analyse, within 
the Bogoliubov model Q , second- harmonic generation of 
excitations in trapped spheroidal atomic condensates: a 
process recently observed in experiment |3j. We find the 
quantal results are significantly different from hydrody- 
namic predictions 5]. In addition to the quantal fre- 
quency shifts, we observe mode hybridization due to cou- 
pling of degenerate Bogoliubov states and the creation 
of coherent-matter dark states. The hybridization pro- 
cess strongly affects second-harmonic conversion. Two 
modes are strongly excited near resonance, and the cou- 
pling strength, and hence conversion rate, to the breather 
mode is half that given by an equivalent hydrodynamic 
estimate. 

Consider then a finite number of atoms, N, each 
of mass m a , occupying a Bose condensate trapped by 
a spheroidal potential such that the angular frequency 
along the polar z— axis, axis of cylindrical symmetry, is 
u>z, and the corresponding radial frequency is lo t . We de- 
fine the trap aspect ratio as A = Lu z /uj r , so that the trap- 
ping potential can be written: Vt ra p = ^m a ljj 2 (r 2 +\ 2 Z 2 ). 
The symmetry means that, in the limit of small ampli- 
tude excitations, the acoustic equation is separable and 
the axial, z, radial r and angular coordinate 6 provide 
good hydrodynamic quantum numbers 0, Q- The az- 
imuthal angular momentum quantum number is denoted 
by m. 

The lowest- frequency excitations of symmetry m = 0, 
are the quadrupole (u>i) and breathing (L02) modes, and 
in the ^-independent hydrodynamic approximation Q: 
w i,2 = 5^(4 + 3 A 2 ± V9A 4 - 16A 2 + 16). The nonlin- 
ear interactions caused by collisions create mode mixing. 



This coupling is most strongly pronounced when tempo- 
ral phase-matching occurs, that is when second-harmonic 
resonance arises: W2 = 2wi. The corresponding values of 
A, the trap aspect ratio, are thus A r = |v77± 5V145 « 
0.683 and 1.952 Q. 

The effect was first observed in an oblate trap by 
Foot and coworkers Q for the following parameters: 
N ~ 20000 atoms of 87 Rb, T < 0.5T C , with ui r = 
2tt x 126Hz. In measurements carried out within the 
range 1.6 < A < 2.8, strongly-enhanced second-harmonic 
generation was found when A = 1.93 ± 0.02. The cor- 
responding chemical potential of this condensate /x, such 
that /x 3> Tiuj r , lie within the range where hydrodynamic 
theory should be valid. And indeed both the frequency 
of the breathing mode, us, Q and the resonant aspect 
ratio, A r , measured by experiment are in very good agree- 
ment with linear hydrodynamic theory 8] . Subsequently 
a hydrodynamic model of the nonlinear mixing of quasi- 
particles was applied to estimate the rate of second- 
harmonic generation and gave results consistent with 
experiment. However, this consistency hides potentially 
important quantal features not previously considered or 
observed. This Letter examines the process of mode mix- 
ing and harmonic generation using a detailed quantal 
treatment of the process Q taking into account atom 
number and quantum pressure corrections. 

The formalism for wave-mixing processes in Bose con- 
densed gases was developed by Burnett and coworkers 3 . 
At low temperatures, T/T c < 0.5, the dynamics of Bose 
condensed gases are dominated by single quasi-particle 
excitations. The atom interactions are represented by 
an s— wave pseudopotential corresponding to a scattering 
length, a s . The spectrum and mode densities of collective 
excitations can be obtained from the Hartree variational 
principle: 



6 J alt (V, [H + yip*ip - ihd t ]ip) = (1) 
where g = (Awh 2 /m a )Na s , H Q = -(h 2 /2m a )V 2 + V tte , p - 
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FIG. 1: Frequency spectrum of the the lowest four even- 
parity, m—0, excitations. The frequencies are in units of the 
radial frequency, u) T - Bogoliubov theory for N = 4650 atoms, 
C = 1000, is compared with hydrodynamic predictions (inset) 
(C —* co). The labelling of the curves in the main graph cor- 
responds to the quantum numbers for a spherical trap: A = 1. 
The frequency of the (+) hybrid mode is marked by the short- 
dashed line, and the (— ) hybrid mode by the long-dashed line. 
The lowest frequency is the I — 2, m — quadrupole (solid 
line). Hybridization of the Bogoliubov breather mode I — 
occurs near A = 0.80, 1.63, 2.60. Second-harmonic resonance 
between the quadrupole and (+) hybrid occurs at A = 1.989. 
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FIG. 2: Hybridization and formation of the Bogoliubov dark 
state. The normalized quasiparticle amplitudes u(r, z) are 
plotted as a function of radial coordinate r (horizontal axis) 
and axial coordinate z (vertical axis) for C = 1000. Red in- 
dicates positive values and blue negative values. Values of 
A increase from left to right: A = 1.35, 1.60, 1.65, 1.75. Be- 
low the crossing, at A=1.35, the surface (upper) and breather 
(lower) modes are distinct. At A = 1.65 the modes combine to 
form the dark state (+) hybrid. As A increases the monopole 
symmetry transfers to the upper(+) hybrid (Fig. 



/i, and (i plays the role of a Lagrange multiplier. The 
condensate and excited modes can be described by the 
linear response ansatz Q where: 

^(r,t) - 6o(*)^(r) (2) 
+ E [&j(*)«i(r)e- fa * t + b*(t)v*(r)e +i ^] 

j>0 

and where <j> represents the highly-occupied condensate; 
that is, |&o| ~ VN > \bj\, j > 0. From the variation 6<j)*, 
and linear expansion in the small parameters bj , b* taken 
as constant, the stationary Gross-Pitaevskii equation and 
Bogoliubov equations follow: 

#o^ + $MV = (3) 

with (4>,4>) = 1. The Bogoliubov modes are solutions of 
the coupled linear equations: 

(Ho + 2g\<j)\ 2 )uj + gcj) 2 v 3 = +huj 3 u 3 (4) 
(H + 2g\(f>\ 2 )v j + g4>* 2 Uj = -hcjjVj (5) 

We take the conventional normalization of modes: 
(ui,Uj) — (vi,Vj) = 5ij. Although the spectrum uij is 
unique, the mode amplitudes may contain arbitrary com- 
ponents of the kernel, the U) — component. After 
diagonalization, it is convenient, though not necessary, 
to enforce orthogonality 0. Gram-Schmidt orthogonal- 
ization ensures the eigenmodes do not overlap the con- 
densate: Uj — Uj — (4>,Uj)4> and Vj = Vj + (<$>,Uj)4>* . 
Suppose that the condensate (mode 0) contains N par- 
ticles, and a single quasiparticle mode i of excitation is 
weakly populated. Then mode-mixing collisions of the 
type, i + i — > + j, and all crossing symmetries, will be 
significant in redistribution of population if phases (en- 
ergies) match and/or the scattering amplitude is large. 

Allowing for variation of the constants bj of the trial 
function © in the variational principle 0J, and neglect- 
ing transitions far from resonance, then we get Q 

ih ^[ = » W <A'<,< ' A ' (6) 

3 

where Ay = uij — 2tOi is the detuning. The mode con- 
version equations are exactly analogous to those in clas- 
sical nonlinear optics, allowing for two or more excited 
modes (j) to be populated from the pump mode (i). The 
coupling strengthjscattering amplitude) for the process, 
Mij, is given by g: 

Mij= 2 fdrW^vtuj+vtvtVj] 

+<j ) [2utvtvj+u* i u* i Uj]} (8) 

and for a condensate with uniform phase throughout, all 
modes can be written in term of real functions. 

The quasiparticle amplitudes and excitation spectra 
were found through discretization of the set of equations 



3 



0.0x10' 



§< -2.0x10"' 
fa 



S -4.0x10"' 



o -6.0x10 : 



o 
O 



3.0x10' 



-10.0x10' 







i i ■ 

i \ 


1 1 l«' 




\ h 












, *t 

,' i 
i 
i 
i 


i| 

ir 

!| 






\ 


i 






\ 

\ 


i 

i, 








\ 1 

v,; 

i j 


i.i. 


1,1, 



0.5 1 1.5 2 2.5 3 

Trap aspect ratio: X=co /co 



FIG. 3: Coupling strength rriij from the m — quadrupole 
state for C = 1000 as a function of trap aspect ratio A: Cou- 
pling to (+) hybrid mode (short-dashed line), to the ( — ) 
hybrid mode (long-dashed line), and to the n = 1,1 = 2 
state (dot-dashed line). The hydrodynamic calculation for 
the breather mode is shown as the dotted line. Second- 
harmonic resonance for the (+) hybrid mode occurs near 
A=0.68 and A=1.95 according to hydrodynamic theory ,8| and 
fig.Q as indicated by the vertical lines. 



(|3I4I5[) using the method of [TO) . We define a dimension- 
less interaction parameter, proportional to the number 
of atoms condensed: C = 8TrNa s (h/2m a uj r )~2, so that 
C — > represents the ideal gas limit. In the hydrody- 
namic approximation, fi h = (15CA/647r) 2 / 5 ?iw r , and the 
hydrodynamic regime can be very roughly characterised 
by x = (/•* /Iv^r) 1. An example of the spectra of Bo- 
goliubov states for C = 1000 ( that is x ~ 7) is shown in 
figEJ This would correspond, for example, to N ~ 4500 
atoms of 87 Rb with scattering length a ~ llOao within 
a trap of frequency uj t = 2-k x 126Hz. Inset, for ref- 
erence, are the equivalent iV-independent hydrodynamic 
results. The main feature of the quantal spectrum, com- 
pared with the hydrodynamic results, are the very large 
differences in the high-Z state frequencies. We find that in 
the quantal regime these surface (high-Z) modes, in par- 
ticular the state labeled I = 4, interact with the second- 
harmonic breather mode and play an important role in 
harmonic generation. 

For low-n, I states, the quantum corrections for the fre- 
quencies are rather small for this number of atoms. For 
example, at A = 1 the solution of equations (|4I5I) for the 
quadrupole state gives uj\ — 1.4808w r for C = 1000. This 
changes to ui% = 1.457w r for C = 2000, and reaches the 
hydrodynamic limit (C — > oo) at lo\ = 1.4142w r . The 
n = 1, 1 = breather frequency is L02 = 2.2101oj r for 
C = 1000, shifting to uj 2 = 2.2195 for C = 2000, tending 
to LO2 = 2.2361clV for C — > 00. However while the quan- 
tal frequencies for these low-lying modes agree well with 
hydrodynamic theory for C ~ 1000, the agreement does 
not extend to wave-mixing process where quantal effects 
are much more pronounced. 

The quantal breather state undergoes a series of 
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FIG. 4: Mode conversion for C = 1000 and 10% pump 
population: (a) A = 2.00: Population of fundamental (n = 
0, 1 = 2, m = 0) quadrupole mode (full line). Population 
of the off-resonant ( — ) mode (n = 0, 1 — 4, m = 0) (long- 
dashed line). Population of resonant second-harmonic (+) 
mode n = 1, 1 — 0, m = (short-dashed line) 
(b) Mode conversion for A = 1.65. The (+) dark-state hy- 
brid has negligible conversion: Population of fundamental 
(quadrupole) mode (full line). Population of (+) dark-state 
hybrid (short-dashed line). Population of (— ) hybrid (long- 
dashed line). 



avoided crossings with higher-Z states as A varies (figure 
nj. The frequency of the I = 4 quantal modes is signif- 
icantly different from the hydrodynamic predictions for 
this regime. For example, at A = 1 the I = 4, m = mode 
has angular frequency lo — 2A26ui r (C — 1000), and 
uj = 2.291w r (C = 2000) compared with u> — 2.000w r as 
C -> 00. While at A = 2, ui = 3.164w r (C = 1000), and 
uj = 3.042w r (C = 2000), compared with uj = 2.732w r 
as C — > 00. As a result, hydrodynamic theory does not 
predict a degeneracy of the / = 4 and monopole states. 

Near the avoided crossings in figure the Bogoliubov 
modes mix symmetries. The hybridization of the modes 
is illustrated by the the excitation functions, u(r, z), 
shown in fig[2 a key finding of this Letter. Below the 
crossing, at A = 1.35, the I = and I = 4 amplitudes 
resemble the spherical pattern and the ^-labelling is cer- 
tainly appropriate. The radial density of the / = 4 mode 
is dominated by the presence of the centrifugal barrier 
pushing the mean radius towards the surface of the con- 
densate. At A = 1.60 (fig. [21 the I — and I = 4 states 
hybridize and lose the character of conventional classifica- 
tion schemes 0,0- Above the degeneracy the modes sep- 
arate and regain their character. Identical features were 
found at the other crossings near A = 0.80 and A = 2.60 
corresponding to interaction with / = 4 and n = 1, 1 = 2 
states, respectively. The hybridization, although primar- 
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ily a quantal feature, persists for much larger numbers 
of atoms. Considering C = 2000, that is N ~ 10 4 , the 
second- harmonic resonance occurs at A r = 1.980 and we 
find a slight displacement of the avoided crossing from 
A c = 1.65 to A = 1.54, but hybridization of the type 
discussed is still present. 

The results for the coupling strengths are shown in 
figure |31 The data are presented in scaled dimension- 
less units in order to compare with hydrodynamic theory 
9J: rritj = 2{h/m a uj r )^ 2 (fi h /huj r )\- 5 ^ M ly The states 
are labelled and identified by the adiabatic noncrossing 
curves in figure^ At the degeneracy A ~ 1.65 the strong 
mixing of states is reflected in the changes in coupling 
strength. Coupling to the upper (+) hybrid mode drops 
suddenly at this point due to destructive interference be- 
tween the mode components. The effect is analogous to 
an optically dark state, that is a coherent superposition 
of states such that the dipole moments cancel [Tij. The 
interference effect is still apparent near A = 2 when the 
state is predominantly breather-like and resonant with 
the quadrupole second harmonic. In contrast the off- 
resonant (— ) hybrid, which is surface-like for A ~ 2 is 
strongly coupled near the degeneracy. 

Similar effects, including the formation of another dark 
state, occur at A w 2.6 corresponding to the crossing be- 
tween I = and I = 6 as shown in figure |3J The result 
of hydrodynamic theory for the quadrupole to breather 
coupling strength [t| is shown in figure |3| as the dotted 
line. The agreement with the quantal calculation is quite 
good near A = 1 but near and beyond the crossing point 
A c = 1.65, the hydrodynamic model does not reflect the 
rapid quantal variation due to hybridization. At the har- 
monic generation resonance near A = 2 the hydrody- 
namic model severely overestimates the coupling strength 
of the breather mode by a factor of two, and neglects the 
contribution of the I — 4 state. 

The solutions of the equations (l(il7l) describing the rel- 
ative populations in modes i and j are well known |12| . 
For a resonant two-state model the characteristic cou- 
pling time is T = gMijb\ (0)|. In terms of the di- 
mensionless matrix elements to? this can be written; 



For A w 1.95, and C = 1000 the quantal prediction is 
|mf 2 | =0.013 compared with the hydrodynamic result 
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=0.028 (the dotted line in figure EJ. Suppose that 
only 5% of atoms were initially in the pump mode, then 
|6J(0)| 2 =0.05. It follows, for C = 1000 and A = 1.95, that 
the respective coupling times are T 9 =14.60 ms compared 
with T h =6.78 ms, which is a substantial and measur- 
able difference. However two excited hybrid modes are 
present near second-harmonic resonance. A numerical 
solution of equations (|6I7I) for a pump population of 10% 
at A = 2 shows the interplay of the hybrid modes: figure 
|H(a). The nonresonant surface hybrid grows faster due 
to the larger coupling strength, but is not as efficiently 
converted as the resonant breather hybrid state. The 
fundamental mode revives Q| after ~ 40ms for these pa- 
rameters. Since the two hybrid modes differ in frequency 
and symmetry this phenomenon should be detectable by 
observation of the radial and axial density variations as 
developed in current experiments |3j. In figure IU(b) cor- 
responding to A = 1.65, very low conversion efficiency is 
observed due to poor phase-matching. Moreover, in this 
case the dark-state hybrid (short-dashed line) is com- 
pletely suppressed. 



In conclusion, we find the hybridization of Bogoliubov 
modes plays a vital role in mode mixing, both in the 
coupling strength and resonant frequency. Another con- 
sequence is the creation of dark coherent states. The 
conversion rate for second-harmonic generation to the 
breather mode is substantially lower than the hydrody- 
namic theory, while the off-resonant surface mode is con- 
verted with almost equal efficiency. This hybridization 
persists for much larger numbers (N ~ 10 4 ) of atoms 
where quantal effects were thought to be neglible. An- 
other potentially interesting case |13| is second-harmonic 
resonant coupling between the to — 2 quadrupole and the 
to — monopole mode. In this case an accidental degen- 
eracy of the monopole with a to = 4 mode occurs exactly 
at the predicted resonance A = ^/16/7 and the effects of 
hybridization would be more extreme. The investigation 
of this particular transition would provide the clearest 
experimental evidence of the effect we have discussed. 
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